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We report density functional studies of the Fei-^Cu^Se alloy done using supercell and coherent 
potential approximation methods. Magnetic behavior was investigated using the disordered local 
moment approach. We find that Cu occurs in a nominal d 10 configuration and is highly disrup- 
tive to the electronic structure of the Fe sheets. This would be consistent with a metal insulator 
transition due to Anderson localization. We further find a strong cross over from a weak moment 
itinerant system to a local moment magnet at x « 0.12. We associate this with the experimentally 
observed jump near this concentration. Our results are consistent with the characterization of this 
concentration dependent jump as a transition to a spin-glass. 
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I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in Fe compounds^ has led to widespread interest both 
in understanding the superconductivity and in unrav- 
eling the normal state properties of these unusual ma- 
terials. These Fe-based materials and cuprates are the 
only known superconductors with critical temperatures 
exceeding 50 K. Therefore it is natural that the relation- 
ship between these two classes of materials has been the 
focus of much investigation. At present, it seems that the 
Fe-based materials are rather different from cuprates at 
least from an electronic point of view - in particular they 
do not show Mott insulating phases, and generally appear 
to manifest more metallic, itinerant electron physics than 
cuprates ^i2il It is unclear whether this means that the 
Fe-superconductors are fundamentally different from the 
cuprates, or whether there is a more subtle connection 
that remains to be found. 

In this regard, the FeSe binary provides a particularly 
useful material for investigation both because of its rela- 
tive simplicity and because of the rich variety of proper- 
ties including superconductivity 5 ^ 6 - found in it, as well as 
in the related FeS, FeTe and alloy systems. FeSe repre- 
sents a robust superconductor with a remarkable increase 
of Tc under pressure reported by many groups ,5i£iL2i2iIfi 
Furthermore, this system is chemically amenable to a 
wide variety of substitutions while still admitting the 
growth of high quality crystalline samples. The response 
of metallic materials to disorder and scattering of various 
types (magnetic, non-magnetic, etc.) induced by substi- 
tutions is a potentially very useful probe of the robustness 
of the metallic phase and of the relationship of physical 
properties and itinerant electron physics. 

The fact that the Fe-based superconductors are less 
strongly correlated than cuprates more readily allows 



the use of standard electronic structure methodology to 
make connections with experiment. However, we note 
that there remain significant errors in the density func- 
tional description of these materials, e.g. in the interplay 
between structure and magnetism, likely due to strong 
spin fluctuations! 11 ' 12 In any case, similar to the other 
Fe-based superconductor s 13 ' 14 ' 15 ' 16 , first principles inves- 
tigations of Fe chalcogenidesii imply very similar physics 
for the pnictide and chalcogenide superconductors. 

Stoichiometry is an important issue for FeSe. The ma- 
terial typically forms with some amount of excess Fe, 
which occupies a crystallographic site outside the Fe 
plane. The formula may therefore be written as Fei+^Se 
or equivalcntly FcSci_j,. Here we use the former nota- 
tion to emphasize the additional Fe site. There have 
been attempts to modulate the excess Fe by alloying with 
Co (by substitution) or Na (by intercalation) ^ Accord- 
ing to that work, the number of excess valence electrons 
(i.e. the experimental doping level of the samples) was 
x ?s 0.027 for Co and x w 0.03 for Na. In both cases 
the Tc was approximately 8.3 K. Other experiments 1 ^ 
showed that x = 0.01 is already sufficient to keep the su- 
perconductivity. At the same time the upper value of x is 
bound by oxygen contamination^ In another work2°- the 
same authors reported that the strongest superconduct- 
ing signal is observed closer to stoichiometry (x = 0.01) 
and that at x = 0.01 the system is much closer to the 
ideal tetragonal phase, while less stoichiometric samples 
have a weaker superconducting transition and magnetic 
contamination. The fact that with x — > the system 
has stronger superconductivity and at the same time be- 
comes magnetically and structurally unstable are indica- 
tions that there may be a quantum critical point (QCP) 
nearby. This view is supported by the fact that near ideal 
stoichiometry Fei.oiSe shows a tetragonal to orthorhom- 
bic distortion as it is cooled through 90Kj21 In other Fe- 
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based superconductors this type of distortion occurs as a 
precursor to the spin density wave (SDW) magnetic or- 
dering (or coincident with the SDW), but Fei.oiSe does 
not show SDW ordering down to the lowest temperature. 

In fact, many authors have discussed the association 
between magnetism and superconductivity in these mate- 
rials, with superconductivity appearing when magnetism 
is suppressed ! 22 ' 23 ! 24 While there is much work pointing 
at the possibility of a QCP affecting the physics of the 
Fe-based superconductors, there is not yet an established 
consensus regarding its nature. On the one hand it has 
been argued that there is a nearness to localization driven 
by Coulomb interactions and that this underlies a quan- 
tum critical point.— However, the phase diagrams do not 
show Mott insulators, and in fact the magnetic phases, 
including the spin-density wave (SDW) phase are unam- 
biguously metallic^ suggesting a quantum critical point 
associated with itinerant magnetism. Experiments prob- 
ing the details of the interplay between magnetic order 
and superconductivity are complicated by the fact that 
chemical disorder in doped samples makes investigation 
of the critical point difficult. As we discuss here, results 
for FeSe indicate that it is very close to the magnetic 
quantum critical point and therefore may be a very use- 
ful system for elucidating the associated physics. 

Before proceeding we mention previous first principles 
studies of defects and doping in iron chalcogenides. Re- 
garding excess Fe, Lee and co-workers early on showed 
that chalcogen vacancies if present will lead to the for- 
mation of ferrimagnetic clusters in FeSe and FeTe^ 
while Zhang and co-workers found that excess Fe in the 
Fei+^Te system serves both as an electron dopant and 
also a magnetic impurity. Han and Savrasov^i reported 
a detailed analysis of the magnetic instabilities and mag- 
netic ordering in relation to the Fermi surface as a func- 
tion of excess Fe concentration. 

The present work is motivated by recent experimen- 
tal findings for Cu-substituted Fei.oiSe^ In particular, 
it was found that Cu substitutes for Fe in the Fe plane, 
and that there is a rapid suppression of Tq in the concen- 
tration range — 4 % and a subsequent metal-insulator 
transition at ~ 4 %. Furthermore, there is a development 
of dynamical magnetic fluctuations detected by NMR, 
which noticeably rise at ~ 12% of Cu substitution. We 
begin with supercell calculations addressing the chem- 
istry and local electronic structure, and then proceed to 
investigate the electronic structure in more detail using 
coherent potential approximation (CPA) calculations for 
the disordered alloy. 



II. SUPERCELL CALCULATIONS 

First principles calculations were performed using 
2x2x1 and 2v / 2x2v / 2xl supercells based on tetragonal 
FeSe (see Fig. [T]) with one Fe replaced by Cu. This 
corresponds to Cu concentrations, x in Fei_a;Cu x Se of 
£=0.125 and x=0.0625, respectively. These calculations 



were done using the generalized gradient approximation 
of Perdew, Burke and Ernzerhof (PBE)>S The lattice pa- 
rameters of the supercells were fixed to the experimental 
values for x=0.12 and £=0.06 as reported in Ref. l28l . 
However, the atomic positions within the cell were de- 
termined by energy minimization. This structural relax- 
ation was performed using the projector augmented wave 
(PAW) method as implemented in the VASP code^ 
with an energy cut-off of 350 cV. The relaxed struc- 
tures showed some expansion of the lattice mainly in the 
in-plane direction around the Cu impurity sites. The 
relaxed Cu-Se bond lengths were ~ 2.44 A, as com- 
pared to Fe-Se bond lengths of 2.25A-2.33A. This is con- 
sistent with the increase in in-plane lattice parameter 
seen experimentally 28 ' 3 1 The electronic structures were 
calculated using the more precise general potential lin- 
earized augmented plane-wave (LAPW) method,-^ with 
the augmented plane wave plus local orbital (APW+LO) 
implementation^ of the WIEN2k co&eM- LAPW sphere 
radii of 2.1 eto for all sites, with well converged basis sets 
and zone samplings were used in this calculation. The 
consistency of the PAW and LAPW calculations was 
checked via the LAPW forces for the structures obtained 
via total energy relaxation with the PAW method. The 
maximum force in the LAPW calculation was 2mRy/ag. 
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FIG. 1: (color online) Crystal structure of Fei-zCu^Se. 

The electronic density of states for the two supercells is 
shown in Fig. while the average Fe and Cu d contribu- 
tions as defined by projections onto the LAPW spheres 
are shown in Fig. [31 As may be seen the Cu d bands are 
at ~ 3 eV binding energy, and are therefore fully occu- 
pied for a nominal d 10 Cu configuration. We also note 
that there are shifts in the Fe density of states as seen in 
Fig. O These shifts are consistent with electron doping. 
However, substitution of a nominal Fc 2+ by Cu 1+ would 
ordinarily be expected to lead to hole and not electron 
doping. This suggests a bonding rearrangement involv- 
ing Se around the Cu site. There is also some Cu d 
contribution to the bands above the Fermi level Ep in- 
dicating covalency involving Cu d states. However, the 
magnitude (sec Fig. [3]) seems inadequate to explain the 
observed electron doping, implying that Se-Se bonding 



3 



6 







1 , t 1 1 

FeSe 

Fe ? CuSe s 


1 ! 


— - Fe, 5 CuSe 16 


! ; 


rf- j '1 


1 1 

ill ! / 




i' 1 A'H 
i 1 KIM 

i ! f in 


" ft \1 1 lr 


1 il i '| j ir'-', n 


i yr mw 


v VZS V,\ 
i ^ A. 

i , 


-6 -4 -2 
E(eV) 


2 



FIG. 2: (color online) Total electronic density of states for 
Fei-zCu^Se supercells on a per formula unit basis. The Fermi 
energy is at eV. 



the arsenides, which leads to the formation of a coherent 
electronic structure with effective doping^ As such, Cu 
alloying cannot be considered with the virtual crystal ap- 
proximation, and more sophisticated treatments such as 
the coherent potential approximation (CPA) 37 ' 38 i 39 are 
needed. Considering the strong scattering implied by this 
result, localization could be caused by disorder, i.e. An- 
derson localization^ should be considered in the context 
of the insulating phase that develops at Cu concentra- 
tions above 4%. Finally, in view of the proximity of these 
phases to magnetism, disruption of the electronic struc- 
ture by scattering (which works against itinerancy) would 
be expected to lead to the formation of local moments 
around the Cu sites. These would be distributed over 
the Fe atoms around the Cu, but because of the d 10 con- 
figuration would not exist on the Cu atoms themselves4i 
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FIG. 3: (color online) Fe and Cu d contributions to the den- 
sity of states for the £=0.125 (top) and £=0.0625 (bottom) 
supercells. These are given by projections onto the LAPW 
spheres for Cu and averaged over the Fe sites, on a per atom 
basis. 



may be important around the Cu site. This is similar 
to findings for CuSe in various structural modifications 
showing an interplay of different bonding types^ 

In any case, from the positions of the Cu d states at 
high binding energy it is clear that Cu substitution is 
highly disruptive to the electronic structure of the Fe 
sheets in FeSe. This is in contrast to the case of Co in 



III. COHERENT POTENTIAL 
APPROXIMATION CALCULATIONS 

Our CPA calculations were performed using the tetrag- 
onal structure (Space group: 115) (see Fig. |T|) , with lat- 
tice parameters from experiment^ similar to the super- 
cell calculations described above. Additionally, in these 
calculations we used experimental atomic positions. It 
follows that with Cu substitution the system linearly 
expands in the plane containing Fe atoms and linearly 
shrinks in the perpendicular direction. The volume in- 
creases weakly with Cu content. At the same time, the 
position of the Se atoms (not fixed by the point group 
symmetry), remains practically constant and is here set 
to z — 0.2707 (in multiples of the c lattice parameter). 
Interestingly, by extrapolating the lattice parameters up 
to 50 % Cu one arrives very close to the structure of 
one of the modifications of CuFeSe2^ That material is 
a semiconductor with low-temperature magnetism i 42 i 43 

The idea of the CPA is to replace the random array 
of real on-site potentials by an ordered array of effective 
potentials. The scattering properties of the effective po- 
tential arc than determined self-consistcntly in terms of 
the local mean-field theory with the condition that the 
total Green's function of the effective system does not 
change upon replacement of the single effective potential 
with the real one. This idea is sketched in Fig. 2] 

SB S3 SB 

FIG. 4: The CPA model: the Green's function of the effective 
medium (gray atoms) is obtained as an average of the partial 
impurities Green's functions (black and white atoms). 



The computational cost of CPA calculations is signif- 
icantly lower than supercell approaches where the disor- 
der is approximated by a randomly generated configu- 
rations within a large number of large supercells. The 
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important advantage is that the CPA does not affect 
the translational (as well as the point group) symme- 
try of the unit cell whereas the supercclls must be chosen 
sufficiently large to avoid effects from the assumed or- 
der or periodic images, which are caused by an artificial 
translational symmetry. The CPA allows investigation 
of the electronic structure as a continuous function of 
the substitution level, which is very important describ- 
ing phase transitions as well as in studying the evolution 
of the electron structure with concentration. The CPA 
also provides type-resolved contributions of the different 
local quantities. The disadvantage is that as any mean- 
field theory, the standard CPA does not include the local 
environment effects such as preferential ordering,— the 
Invar effec t 45 ' 46 and lattice relaxations around the impu- 
rity site. The low concentrations of interest here, and 
the results of the superccll calculations, which find mod- 
est changes in local structure, indicate that the effects of 
local ordering and lattice distortion are not likely to be 
the crucial factors in determining the electronic structure 
of the alloy. 

The Fei_a,Cu x Se system is known to be non-magnetic, 
in the sense of not having ordered magnetism, over the 
investigated range of Cu concentrations. However, as 
mentioned, there is experimental evidence that magnetic 
fluctuations noticeably strengthen above 12% Cu. The 
consistent way of studying these fluctuations is to deter- 
mine the coefficients Ai of the Ginzburg-Landau energy 
expansion E = Epm + X^i AiM 21 in the local magnetiza- 
tion Mj 4 ^ The best local approximation to the disordered 
paramagnetic state £pm is given by the so-called disor- 
dered local moments (DLM) approach ) 48 ' 49 which is used 
in the present work. The DLM effective medium gives 
an accurate representation of a paramagnetic state with 
randomly oriented spins by using an cquiatomic random 
alloy of spin-up and spin-down atoms, in complete anal- 
ogy to the CPA (see Fig. Op. 




FIG. 5: The DLM model: the Green's function of the effective 
non-magnetic medium is obtained as an average of the partial 
polarized impurity Green's functions. 

In this context the DLM picture provides a common 
basis to study different types of possible magnetic or- 
der. The important feature of the DLM is the consis- 
tent description of the paramagnetic state for the systems 
with itinerant and local moment behavior j 47 ' 50 local mo- 
ment systems keep their finite moments in the paramag- 
netic state whereas itinerant systems become locally non- 
magnetic. The essential assumption of the DLM is the 
quasi-static treatment of fluctuations, ft assumes that 
the fluctuations of the local magnetic moment are much 
slower than the electron motion, so that the instanta- 
neous band structure corresponds to the equilibrium elec- 



tronic state. This should lead to a certain overestimation 
in case of large fluctuating magnetic moments and to a 
certain suppression in case of small moments, by sharply 
emphasizing magnetic transitions. 

Here we use the KKR (Korringa-Kohn-Rostoker) 
method in the Munich SPR-KKR (spin-polarized rela- 
tivistic) package^ 1 - for our CPA calculations. These cal- 
culations were done with the local spin-density approxi- 
mation (LSDA) in the form of Vosko, Wilk and Nusair— . 
Here we use the full-symmetry potential (FP) method to 
account for the non-spherical contributions to the one- 
electron potential. This is important because the layered 
structure of FeSe is strongly anisotropic. As expected, 
relativistic effects are found to be relatively small. Thus 
although the calculations are fully relativistic, only spin 
moment is reported in the following as the orbital mo- 
ments are found to be negligibly small. 

The effective DLM medium is represented by the ran- 
dom Feo.sFeo.s alloy, with ± 2 hb anti-parallel spins ini- 
tially induced for each on-site component. The system 
is then allowed to relax during the self-consistent itera- 
tion process. In case of locally non-magnetic solutions 
the proximity of the system to the magnetic state is es- 
timated by considering the Fe local spin susceptibility 
X = dfi I dB e xk | H — o • This is calculated by applying the 
small external local magnetic field of <iB 0Xt = 2 mRy. 

The symmetry of the magnetic fluctuations are con- 
sidered to be characterized by the exchange coupling 
constants Jy of the Heisenberg local moments picture: 
H e s = — X^i>j Jij fiifij ) where pn is the local magnetic 
moment of the i-th atomic site. However, it should be 
emphasized that the use of a Heisenberg model in analyz- 
ing low energy magnetic fluctuations does not necessar- 
ily imply local moment magnetism. Here, the so-called 
real-space formalism is utilized^. Again, in case of the 
non-magnetic solution, a magnetic moment is induced by 
a small external field of 2 mRy. 



IV. CPA RESULTS 
Magnetic moments and susceptibility 

Regardless of the calculation mode (magnetic or non- 
magnetic, fully- or non-rclativistic, spherical or non- 
spherical potential), we find that all properties (den- 
sity of states (DOS) at the Fermi energy, magnetic mo- 
ments, susceptibility, etc.) calculated as a function of x 
exhibit discontinuities or jumps at a Cu concentration, 
x m « 12%. Importantly, this concentration corresponds 
remarkably well to the experimentally observed onset of 
the dynamical magnetic fluctuations^ The DLM cal- 
culations lead to the appearance of local moments at 
the Fe site that vary weakly within a given regime: 
low Cu-concentration (x < x m ) or high Cu-concentration 
(x>x m ) (Fig.EKa)). 

Within the low Cu-concentration regime the magnetic 
moment is suppressed almost to zero (~ 10 _3 /kb), how- 
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Cti concentration x [%] Cu concentration x [%] Cn concentration x \%] 



FIG. 6: Local magnetic moments of Fe (a), local spin sus- 
ceptibility (b) and the total DOS at the Fermi energy (c) 
calculated as a function of Cu concentration. Filled squares 
mark DLM, hollow squares - non-magnetic calculations. 

ever at x = x m it rises sharply to about 2.35 /xb- The 
Fe local spin susceptibility shows a similar discontinu- 
ous behavior at shown in Fig. [HJb). On the 
other hand, besides the sharp magnetic transition at x m , 
the magnetic susceptibility indicates that the system ap- 
proaches magnetism in the low Cu-concentration regime 
as well, however its ground-state remains non-magnetic. 
As may be seen (Fig. [6jc)), the changes in magnetic be- 
havior as a function of x closely follow the behavior of 
the density of states at the Fermi energy and is therefore 
related to the band structure. 

We emphasize that the DLM calculations provide in- 
formation about moment formation in the disordered 
(paramagnetic) state, not directly about magnetic order- 
ing. In fact, at the LSD A or GGA levels FeSe is already 
unstable against magnetic ordering, specifically against 
SDW order ^ This ordering is presumably suppressed by 
spin-fluctuations the detailed nature of which is not fully 
understood at present, but which might be central in un- 
derstanding the superconductivity of these phases. What 
a cross-over from a non-magnetic to a state with moments 
in the DLM indicates is a cross-over from a state that may 
have itinerant magnetism to a state with stable moments 
independent of ordering. This is therefore a change in the 
nature of the magnetism at this value of x. While such a 
cross-over is not directly a QCP it may be an indicator. 
FeSe, as mentioned, does not show long range magnetic 
order, but magnetic ordering is very difficult to avoid in 
clean metals with stable moments and electronic struc- 
tures that arc neither low-dimensional nor geometrically 
frustrated. The effect of the disorder induced by rela- 
tively low amounts of Cu alloying appears to be a change 
of magnetic character from a softer moment (itinerant) 
regime to a regime with more robust local moments. 




Cu concentration x [%] 



FIG. 7: (Color online) Exchange constants J i (i — 1,2,3,4) 
for the 4 nearest-neighbor in-plane couplings as a function 
of Cu concentration. Each Joi coupling curve is marked by 
corresponding index i. 

Since in the range of x < x m , Joi has larger am- 
plitude than Jq2, while the number of the first- and 
second-nearest neighbors is the same, the nearest anti- 
ferromagnetic order seems to be more preferable. For 
the pure FeSe the particular order was already studied 
in more details by other authorsi 17 i 27 ' 54 i 55 In the extrap- 
olation from the low Cu-concentration limit the present 
exchange coupling constants agree well with those pub- 
lished bcforei 17 ' 54 On the other hand, as we have seen, the 
system does not have stable magnetic moments within 

X X J^Q . 

The monotonously increasing amplitude of the nearest- 
neighbor coupling indicates the growing electron dereal- 
ization with increasing Cu content. However, at x — x m 
one observes a rapid suppression of all Jij , which in turn 
indicates the sudden localization of the d-electrons and 
corresponds to the development of the strong atomic mo- 
ments. In this regime Joi ~ Jo2, while the higher-order 
coupling constants nearly vanish. This situation with sta- 
ble local moments emerging at x m 0.12, competing ex- 
change interactions, and disordered Cu positions (which 
correspond to sites with no moment) is consistent with 
the formation of a spin-glass. 



V. BAND STRUCTURE ANALYSIS 



Exchange coupling 

Turning to the exchange couplings, we find very small 
inter-plane interactions below the precision of the cal- 
culations. Furthermore, the in-plane contributions are 
significant only for the first four in-plane coordination 
circles. The corresponding values are shown as a func- 
tion of Cu concentration in Fig. [7J 



Band structure 

In the following we consider the Bloch spectral func- 
tions, which take the role of the band structure in the dis- 
ordered system. These are plotted along the path shown 
in Fig. [8l The calculated DOS and the corresponding 
spectral functions arc shown in Fig.[9J The general trend 
in the low Cu-concentration regime is the redistribution 
of the spectral weight from the M- A onto the Z-T line, by 
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shifting the whole picture down in energy with respect to 
the Fermi level. Qualitatively, this is exactly what was 
found in our supercell calculations where in spite of the 
Cu d 10 configuration, Cu substitution produces electron 
doping. 
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coordinates are given in units of (2n/a, 2n/a, 2n/c). 
marked details are discussed in the text. 
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As in the other Fe-based superconductors the nearness 
to magnetism of FeSe is related to the electronic struc- 
ture specifically the high density of states derived from 
Fe d orbitals and the nested Fermi surface. The Fermi 
energy at x=0 occurs towards the bottom of a pseudo- 
gap in the DOS. As seen in the Bloch spectral functions 
there are two effects as Cu is alloyed into the material. 
First of all the Fermi energy is raised and secondly, as 
a consequence of the disorder, the DOS as a function of 
energy is smoothed, filling in the pseudogap. This ex- 
plains the initial decrease, followed by an increase in the 
tendency towards moment formation. When the Cu con- 
centration reaches x m value there is a dramatic change 
in the spectral function. Since the system develops large 
local moments, the DLM spin disorder leads to a strong 
incoherence in the bands. 



Fermi surface 

To study the nesting effects we track the substitution- 
induced transformation of the Fermi surface in the BZ. 
In the following we will consider the Fermi surface cross- 
sections by horizontal (k z = const) and vertical (k x = k y ) 
planes. 

In general, the electronic structures of the Fe-based 
superconductors show disconnected Fermi surfaces con- 
sisting of hole sections around the zone center (along T-Z) 
and two electron sections at the zone corner (along M-A 
direction). Our calculated Fermi surface cross-sections 
(Fig. UOp arc in good agreement with this picture. 

It turns out that the growing chemical disorder with 
increasing Cu-concentration destroys the Fermi surface, 
especially affecting the hole-pockets (centered at T-Z). 
Since the in-plane hole-electron nesting is the main mech- 
anism for the superconductivity in FeSe and related 
compounds^ the superconducting signal must strongly 
attenuate. Indeed, the experimentally measured Tq falls 
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down very fast and completely vanishes at x 4 % 
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FIG. 9: (Color online) Atomic type-resolved DOS and the 
corresponding Bloch spectral functions for x — 1, 11, 12 % as 
calculated within the DLM. 



Although the reason for the magnetic transition can 
be understood in terms of the band structure (Fig. [9]), 
it is still useful to observe the accompanying changes of 
the non-magnetic Fermi surface. The corresponding re- 
sults for 11 % and forced non-magnetic calculations for 12 
and 14 % are compared on Fig. [TTJ Here the low-weight 
contributions are cut off, otherwise the band structure 
analysis will be much complicated by chemical disorder. 

Thus, the Fermi weight rapidly "spills" from k z = 0.5 
onto the k z = plane. This shows a loss of the Fermi sur- 
face in the k z = 0.5 plane first as a function of Cu alloying 
although there is a strong loss of spectral weight at all k z 
consistent with insulating character. At the same time 
the central figure does not blow up any more, but rather 
pulls the borders of the magnetic BZ (MBZ) (the two rel- 
evant in-plane MBZs are marked by the dashed and the 
dotted lines on Fig. ITT]) ■ Analogous behavior was shown 
for the related Fei +x Te compound^: it was suggested 
that such redistribution of the Fermi weight between two 
neighboring MBZs reflects the competition between mag- 
netic interactions of different symmetries. The inhomo- 
geneity becomes especially sharp at x — 14 %, indicating 
that the system starts to prefer the certain in-plane mag- 
netic order for x > x m . 

Preferential accumulation of the Fermi weight on the 
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FIG. 10: Cross-sections of the Fermi surface in BZ by k z — 0.5 
(a, d, g), k x = k y (b, e, h) and k z = (c, f, i) planes. 
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FIG. 11: Cross-sections of the Fermi surface in BZ by k z — 0.5 
(a, d, g), k x = k y (b, e, h) and k z — (c, f, i) planes. Calcu- 
lations are done within the non-magnetic regime. 



edges of the corresponding MBZ in the k z = plane indi- 
cates the derealization of the conduction electrons in real 
space. This derealization becomes unfavorable, since 



the in-plane lattice spacing increases (the accompanying 
vertical real-space derealization is a minor effect). At a 
certain point (x = x m ) the system undergoes the localiza- 
tion transition. As a result, the strong magnetic moment 
arises. This rapid localization is also reflected by the re- 
duced amplitude of the exchange coupling constants for 
x > x m as discussed above. 



VI. SUMMARY AND CONCLUSIONS 

To summarize, we find that Cu occurs in a d 10 con- 
figuration when alloying into FeSe. Nonetheless it serves 
as an effective electron dopant. Importantly, it is also 
a source of strong scattering. The calculated changes in 
the Bloch spectral function with alloying show no signs 
of gapping near the Fermi level over the concentration 
range studies. This situation suggests that the insulat- 
ing phase occurring above ~ 4% Cu may be an Anderson 
localized system arising from disorder rather than a con- 
ventional semiconductor. The magnetic instability ob- 
served at about 12% Cu is characterized as a spin-state 
like transition, where a soft moment near magnetic state 
at low x gives way to a state with stable disordered local 
moments at high x. This is consistent with the formation 
of a spin-glass especially considering the near compensa- 
tion of nearest and next nearest exchange interactions. 
This is in accord with experimental observations^. 

The result that the Fermi weight for the electron sec- 
tions at k z = 0.5 is suppressed before that at k = 
suggests a three dimensional character to the electronic 
structure and scattering above x ~ 0.1. It will be quite 
interesting to investigate this possibility experimentally, 
via single crystal transport or experiments under strain, 
especially uniaxial strain if this becomes feasible. 

Finally, we note that the present calculations assume 
full disorder between Cu and Fe, and also that both Fe 
and Cu remain in the Fe plane and do not occupy other 
sites, such as the excess Fe site in the interstitial. It will 
be useful to test the extent to which these assumptions 
hold for the real material, for example by using diffraction 
to refine the occupancy of the interstitial site. 



Acknowledgments 

The authors are much indebted to Roser Valenti and 
Jurgcn Kiibler for helpful discussions. The financial 
support by the SFB/TRR49 is gratefully acknowledged. 
Work at Oak Ridge National Laboratory was supported 
by the Department of Energy, Division of Materials Sci- 
ences and Engineering. 



* Electronic address: chadov@uni-mainz.de 



Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 



8 



Amer. Chem. Soc. 130, 3296 (2008). 

2 S. E. Sebastian, J. Gillett, N. Harrison, P. H. C. Lau, D. J. 
Singh, C. H. Mielke, and G. G. Lonzarich, J. Phys. Con- 
dens. Matter 20, 422203 (2008). 

3 D. H. Lu, M. Yi, S. K. Mo, A. S. Erickson, J. Analytis, 
J. H. Chu, D. J. Singh, Z. Hussain, T. H. Geballe, I. R. 
Fisher, et al., Nature (London) 455, 81 (2008). 

4 W. L. Yang, A. P. Sorini, C. C. Chen, B. Moritz, W. S. Lee, 
F. Vernay, P. Olalde-Velasco, J. D. Denlinger, B. Delley, 
J. H. Chu, et al., Phys. Rev. B 80, 014508 (2009). 

5 F.-C. Hsu, J.-Y. Luo, K.-W. Yeh, T.-K. Chen, T.-W. 
Huang, P. M. Wu, Y.-C. Lee, Y.-L. Huang, Y.-Y. Chu, 
D.-C. Yan, et al., Proc. Nat. Acad. Sci. (USA) 105, 14262 

(2008) . 

6 Y. Mizuguchi, F. Tomioka, S. Tsuda, T. Yamaguchi, and 
Y. Takano, Appl. Physics Lett. 93, 152505 (2008). 

7 S. Medvedev, T. M. McQueen, I. Troyan, T. Palasyuk, 
M. Eremets, R. J. Cava, S. Naghavi, F. Casper, V. Kseno- 
fontov, G. Wortmann, et al., Nature Materials 8, 630 

(2009) . 

8 S. Margadonna, Y. Takabayashi, Y. Ohishi, Y. Mizuguchi, 
Y. Takano, T. Kagayama, T. Nakagawa, M. Takata, and 
K. Passides, Phys. Rev. B 80, 064506 (2009). 

9 G. Garbarino, A. Sow, P. Lejay, A. Sulpice, P. Toulemonde, 
W. Crichton, M. Mezouar, and M. Nunez-Regueiro, Euro- 
phys. Lett. 86, 27001 (2009). 

10 T. Imai, K. Ahilan, F. L. Ning, T. M. McQueen, and R. J. 
Cava, Phys. Rev. Letters 102, 17705 (2009). 

11 I. I. Mazin, M. D. Johannes, L. Boeri, K. Koepernik, and 
D. J. Singh, Phys. Rev. B 78, 085104 (2008). 

12 D. J. Singh, Physica C 469, 418 (2009). 

13 L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. Rev. 
Letters 101, 026403 (2008). 

14 Z. P. Yin, S. Lebegue, M. J. Han, B. P. Neal, S. Y. 
Savrasov, and W. E. Pickett, Phys. Rev. Letters 101, 
047001 (2008). 

15 D. J. Singh and M. H. Du, Phys. Rev. Lett. 100, 237003 
(2008). 

16 I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 

17 A. Subedi, L. Zhang, D. J. Singh, and M. H. Du, Phys. 
Rev. B 78, 134514 (2008). 

18 Z. Liu, A. Fang, F. Huang, and M. Jiang, 
arXiv:0808.1784vl (2008). 

19 T. M. McQueen, Q. Huang, V. Ksenofontov, C. Felser, 
Q. Xu, H. Zandbergen, Y. S. Hor, J. Allred, A. J. Williams, 
D. Qu, et al., Phys. Rev. B 79, 014522 (2009). 

20 A. J. Williams, T. M. McQueen, and R. J. Cava, Solid 
State Commun. 149, 1507 (2009). 

21 T. M. McQueen, A. J. Williams, P. W. Stephens, J. Tao, 
Y. Zhu, V. Ksenofontov, F. Casper, C. Felser, and R. J. 
Cava, Phys. Rev. Lett. 103, 057002 (2009). 

22 C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. R. II, J. L. 
Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. L. Wang, 
et al., Nature (London) 453, 899 (2008). 

23 G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, 
P. Zheng, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 
100, 247002 (2008). 

24 F. Ma and Z. Lu, Phys. Rev. B 78, 033111 (2008). 

25 J. Dai, Q. Si, J.-X. Zhu, and E. Abrahams, Proc. Nat. 
Acad. Sci. USA 106, 4118 (2009). 

26 K.-W. Lee, V. Pardo, and W. E. Pickett, Phys. Rev. B 78, 
174502 (2008). 

27 M. J. Han and S. Y. Savrasov, Phys. Rev. Letters 103, 



067001 (2009). 

28 A. J. Williams, T. M. McQueen, V. Ksenofontov, C. Felser, 
and R. J. Cava, J. Phys.: Condens. Matter 21, 305701 
(2009). 

29 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

30 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

31 Y. Li, X. Lin, Q. Tao, C. Wang, T. Zhou, L. Li, Q. Wang, 
M. He, G. Cao, and Z. Xu, New J. Phys. 11, 053008 (2009). 

32 D. J. Singh and L. Nordstrom, Plane-waves, Pseudopoten- 
tials and the LAPW Method, 2nd Ed. (Springer Verlag, 
2006). 

33 E. Sjostedt, L. Nordstrom, and D. J. Singh, Solid State 
Commun. 114, 15 (2000). 

34 P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2k, An Augmented Plane Wave + Local 
Orbitals Program for Calculating Crystal Properties (K. 
Schwarz, Tech. Univ. Wien, Austria) (2001). 

35 V. Milman, Acta Cryst. B58, 437 (2002). 

36 A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, D. J. Singh, 
and D. Mandrus, Phys. Rev. Lett. 101, 117004 (2008). 

37 B. Gyorffy, Phys. Rev. B 5, 2382 (1972). 

38 G. M. Stocks and H. Winter, The electronic structure of 
complex systems (Plenum Press, New York, 1984), p. 463. 

39 W. H. Butler, Phys. Rev. B 31, 3260 (1985). 

40 P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

41 Another interesting possibility is that moments could form 
on the neighboring Fe but that these could form singlets 
where the four Fe atoms around the Cu could have can- 
celling moments due to antiferromagnetic interactions be- 
tween them. Recalling that the four Fe atoms around a 
Cu substitutional are second neighbors on the Fe square 
lattice, the case where the Fe moments are aligned would 
correspond locally to the checkerboard (nearest neighbor 
antiferromagnetic) order, while the antiferromagnetic state 
with cancelling Fe moments would correspond locally to 
the SDW pattern. 

42 Y. A. Dorofeev, A. Z. Men'shikov, and V. G. Pleshchev, 
Physics of the Solid State 36, 1329 (1994). 

43 J. C. Woolley, A. M. Lamarche, G. Lamarche, R. Brun del 
Re, M. Quintero, F. Gonzales- Jimenez, I. P. Swainson, and 
T. M. Holden, J. Magn. Magn. Mater. 164, 154 (1996). 

44 S. Mankovskyy and H. Ebert, Phys. Rev. B 74, 54414 
(2006). 

45 F. Baudelet, S. Odin, C. Giorgetti, E. Dartyge, J. P. Itic, 
A. Polian, S. Pizzini, A. Fontaine, and J. P. Kappler, J. 
Phys. (Paris) 7, C2 441 (1997). 

46 V. Crisan, P. Entel, H. Ebert, A. Akai, D. D. Johnson, and 
J. B. Staunton, Phys. Rev. B 66, 014416 (2002). 

47 N. M. Rosengaard and B. Johansson, Phys. Rev. B 55, 
14975 (1997). 

48 A. J. Pindor, J. Staunton, G. M. Stocks, and H. Winter, 
J. Phys. F: Met. Phys. 13, 979 (1983). 

49 B. L. Gyorffy, A. J. Pindor, J. Staunton, G. M. Stocks, 
and H. Winter, J. Phys. F: Met. Phys. 15, 1337 (1985). 

50 A. V. Ruban, S. Khmelevskyi, P. Mohn, and B. Johansson, 
Phys. Rev. B 75, 054402 (2007). 

51 H. Ebert and M. Battocletti, Solid State Commun. 98, 785 
(1996). 

52 S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 
1200 (1980). 

53 A. I. Liechtenstein, M. I. Katsnelson, V. P. Antropov, and 
V. A. Gubanov, J. Magn. Magn. Materials 67, 65 (1987). 



54 Y.-F. Li, L.-F. Zhu, S.-D. Guo, Y.-C. Xu, and B.-G. Liu, Letters 102, 177003 (2009). 

J. Phys.: Condens. Matter 21, 115701 (2009). 



or, 



F. Ma, W. Ji, J. Hu, Z.-Y. Lu, and T. Xiang, Phys. Rev. 



